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Abstract 

The light-front coupled-cluster method is a technique for solving Hamiltonian eigenvalue prob- 
lems in light-front-quantized field theories. Its primary purpose is to provide a systematic sequence 
{3jQ[ of solvable approximations to the original eigenvalue problem without the truncation of Fock space. 

Here we discuss the incorporation of zero modes, modes of zero longitudinal momentum, into the 
formalism of the method. The use and importance of these modes are illustrated in two-dimensional 
model field theories. 
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I. INTRODUCTION 



One approach to the nonperturbative solution of a strongly interacting quantum field 
theory is that of Hamiltonian methods in light-front quantization [j], 0]. The Hamiltonian 
eigenvalue problem is formulated in Fock space as the fundamental equation 

/If 2 -l- P 2 

v-mp))= p + m ^(p)), (i.i) 

where V~ is the light-front energy operator, conjugate to the light-front time x + = t + z, 
and P_ = (P + = E + P z , P± = (P x , P y )) is the light-front momentum, with P + conjugate to 
x~ = t — z. The eigenstate has mass M and is expanded in a Fock basis of eigenstates 
of P and of particle number. The coefficients of the Fock states are the wave functions that 
describe the eigenstate. The eigenvalue problem (11 .1ft is equivalent to an infinite coupled 
system of integral equations for these wave functions. 

Light-front quantization is quite convenient for this construction. In particular, the Fock- 
state expansion is well-defined, and there is a separation of external and internal momenta 
for the constituents, which leads to well-defined and boost-invariant wave functions It 
does, however, have the disadvantage that, except for modes of zero longitudinal momentum 
P + |3H5|, there are no contributions to the vacuum, making the perturbative vacuum JO) the 
physical vacuum. The calculation of vacuum effects, such as symmetry breaking jj-l9|, is 
then a challenging aspect of light-front quantization, and the incorporation of zero modes 
can be important in some calculations. Here we consider zero modes in the context of the 
light-front coupled-cluster (LFCC) method 11 



To have a finite calculation for an eigenstate, the usual step is a truncation of Fock space, 
to have a finite number of wave functions and a finite set of equations. This, however, 



leads to many difficulties, particularly uncanceled divergences [12J. The LFCC method is 
designed to avoid these difficulties by not truncating Fock space but instead restricting the 
relationships between wave functions in such a way as to produce a finite set of (nonlinear) 
equations. The mathematical structure of the method is closely related to that of the many- 



body coupled-cluster method [L3( used in nuclear physics and physical chemistry \U 

The LFCC method constructs the eigenstate \ip) from a valence state \<p) with the smallest 
number of constituents and the exponentiation of an operator T that increases the particle 
numbers, to generate higher Fock states. The general form is 

|V(P)) = v^e T |0(P)), (1.2) 



with yZ a normalization factor. The eigenvalue problem is then converted to a valence 
eigenvalue problem 

\l- _|_ P2 

p v v-\m) = p + - 1 10(g)), (i-3) 

with V~ = e~ T V~e T an effective Hamiltonian and P v a projection onto the valence sector, 
and to an auxiliary equation for T 

(1- J R U )^|0(P)) = O. (1.4) 

The auxiliary equation is actually an infinite set of equations for the infinite set of terms in 
T, and as such we still have an exact representation of the original eigenvalue problem. The 
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approximations that lead to a finite set of equations, without truncating Fock space, are 
the truncation of T to a finite number of terms and a matching truncation of the projection 
1 — P v to generate the finite number of equations needed to solve for the terms in T. 

Here we wish to take a step toward the inclusion of vacuum effects, such as spontaneous 
symmetry breaking in wrong-sign 4 , by introducing zero-mode contributions to the T oper- 
ator in the LFCC method. We do this as a limiting procedure, with modes of momentum e 
introduced at the start of a calculation and the limit e — > taken at the end. The technique 
is developed in a series of examples; we discuss 3 theory [l5j], 4 theory [l6l- 18|. and the 



Wick-Cutkosky model |19| in Sees. II, III, and IV, respectively. A brief summary is given 
in Sec. |V] Some details for singular terms are included in an Appendix. 



II. 3 THEORY 

The Lagrangian of 3 theory is 

C= 1 -(d, ( pY~ 1 -^-^. (2.1) 
From this, the two-dimensional light-front Hamiltonian density is 

H = \d^d + <P -£ = + t^ 3 - (2-2) 

The mode expansion for the field at zero light-front time is 

dp+ ( a(p + )e- ip+x ~/ 2 + a\p + )e ip+x ~' 2 ) , (2.3) 



with the modes quantized such that 

[a(p + ), a\p' + )} = 5(p + - p' + ). (2.4) 
The normal-ordered light-front Hamiltonian V~ = V^ ee + V^ t is then specified by 

Kee = J dp + ^(p + )a(p+) (2.5) 



and 



P " = ~2 I T^faj / [«VW^+W0 - (2-6) 

+a) (xp + )a i ((1 — x)p + )a(p + )~\ . 



The form of V^ t , in terms of the momentum p + flowing through the operator and the 
longitudinal momentum fraction x, is particularly convenient for what is to follow. For 
a graphical representation of V~ , see Fig. [TJ The effective Hamiltonian V~ of the LFCC 
method is computed from the Baker-Hausdorff expansion 

= e- T V~e T = V~ + [P-,T} + -[[P-,T],T} + --- , (2.7) 



3 



- -x- - > »C 




(b) 

FIG. 1. Diagrammatic representation of (a) the light-front Hamiltonian V~ and (b) the approx- 
imate T operator for 3 theory. A circle with a v represents a zero-mode contribution; a circle 
with a 1 represents the function t\ found in the T operator. The cross represents the kinetic en- 
ergy contribution. External lines on the right represent annihilation operators; those on the left, 
creation operators; and internal lines, contractions. 



given an approximation for T. 

For <f) 3 theory, the simplest approximation for T is a single boson emission process repre- 
sented by 

Ti = J dx J dp + \j^-t l (x)a\xp + )a\{l-x)p + )a{p + ). (2.8) 

Here ti(x) is a function to be determined by the auxiliary equation (jl.4p . From the structure 
of T we can see that only the symmetric part of t\ contributes; therefore, we can assume 
t\(x) — — x). The corresponding truncation of the projection 1 — P v is to include only 
states with one particle in addition to those in the valence state. 

To include the zero mode, we introduce another operator to T = T + T\ as 

T = A e a a^(e). (2.9) 

This creates a particle of momentum e, which does violate momentum conservation for T; 
however, the limit of e — > will be taken at the end, so that To is associated only with zero 
modes. The power a will be determined from the auxiliary equation, but the outcome can 
be anticipated by considering the corresponding transformation of the field: 

e- T °0 e ^o = + [0, r ] = + 4=, (2.10) 



which indicates a = 1/2 and a shift of v = A / yAw for the field. A graphical representation 
of T is given in Fig. [TJ 

A calculation of the terms in the Baker-Hausdorff expansion (12. 7ft then determines the 
effective Hamiltonian. Only a finite number of terms will contribute to the eigenvalue 
problem, because we need only terms that change the particle number by no more than one; 
any more than this would go beyond the truncation of the projection (1 — P v ). We will 
consider only the vacuum and the one-particle state as valence states, and, therefore, terms 
in V~ with more than one annihilation operator will also be neglected; however, terms with 
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more than one do need to be kept in intermediate calculations of commutators, because 
higher-order commutators can reduce the total number of annihilation operators. 
We compute the following commutators for To: 

e 

[[P^,T ],T ] = 0, 



[K nV T ] = X J 



dp^ 



A e a 



v/4ttp+ ^/ e (p+ - e) 



w (p + )a(p + — e), 



(2.11) 
(2.12) 

(2.13) 



[[Vr v T ],T } = \ 



a f (e). 



For [V int ,T ], we keep only the leading term as e — > 0; this leaves 

[Pr t ,T ]^x[ %^=a\p+)a{p+\ 
J P^ yAite 



(2.14) 



(2.15) 



which we use in the computation of [[P^ v T ], T ]. There is a subtle aspect of the e — > 
limit, in that the integration range of p + includes a region where it is not large compared to 
e; a more careful calculation is described in the Appendix, with the same result. 
The commutators for 7\ are 



dp A 



[V int ,T 1 } = J dp + -La\p + )a{p + 

dx 



x(l — x) 
x a\xp + )a) '((1 - x)p + )a(p + ), 



(2.16) 



(2.17) 



\/x(l — x 

X 



dx'ti(x') 



^a\x'p + )a\(l - x')p + )a(xp + )a((l - x)p + ) 

4:71 

1 — x f dp^ 



+2A / ; dx [ dx'tAx')^^-^- [ 

J y/x(l- X )J U ll-X'J 

dp'+ 



47T 



x 



5{{l - x)p + - (1 - x')p' + )a\p + )a\x'p' + )a(xp + )a(p' + ), 



[[V^,T 1 ],T 1 ]= dx 



dp A 



^4ttp+ 

with the self-energy integral 



h(x) 



-21 



+ V(x) 



x(l — x) 

x a\xp + )a\{l — x)p + )a(p + ), 



(2.18) 



I = — 



dxti(x) 



and vertex function 



V(x) 
V(x) 



1 

2 
A 



47rJ y/ x (l - x) 

V(x) + V(l-x) 

ti(x/y)h(y) 



IT 



dy 



\/y{ l -y){y - x){i- x) 



(2.19) 

(2.20) 
(2.21) 
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Here and below we have taken advantage of the symmetry t\{x) = ti(l — x) described 
following Eq. (12.81) . The double commutators that mix Tq and T\ are 



XA n e a 



ire 



dx 



dp^ 



- : ti(x) 



x(l — x) 
x a} (xp + )a) ((1 — x)p + )a(p + ) 



-A e a t 

a 1 

e 



a' e 



(2.22) 



(2.23) 



V4 



7Tf 



-Mix) 



- 1 



x a^(a;p + )a^((l — x)p + )a(p + ). 



That the self-energy integral is the same in all terms is one of the hallmarks of the LFCC 
method. For calculations with Fock-space truncation, the self-energy contributions are sector 
and spectator dependent (l2| . 

From these commutators we build the expression for V~ as 



v- 



A)— 

e 



" + 2 /+ 27S 



aUe) 



(2.24) 



dp A 
p+ 

+ I dx 



a\p + )a(p + ) 



' a^(xp + )a^((l — x)p + )a(p^ 



/i 2 + / + 



AA e c 



A/47re 
1 



x(l — x) 



- 1 



*i(z) 



A 1 Kn ^ 1 

H ; + -y(ar) ; 

2 2 2x l-i 



We have kept only those terms that do not annihilate the one-particle Fock state. A graphical 
representation is given in Fig. [2] 



For a one-particle valence state \<j>) = aJ(P + )\Q), the eigenvalue problem P V V~\4>) = ¥^ 



is 



/i 2 + / + 



7l€ 



M 



p+ IV/ 

2.25) 



The auxiliary equation (jl.4p in general includes two and more particle sectors. For our 
approximation, 1 — P v is truncated to two-particle states only, and (jl.4p becomes 



+ 



^ + n 1 + 



1 AA e a 



2 v/4 



7TC 



a T (e)a t (P + )|0) 



(2.26) 



(it; 



=a t (xP + )a t ((l -x)P + )\0) 



[i 2 + 1 + 



XA^ 



1 



x(l — x) 



t x (x) 
I 



2x(l - x, 



h{x) 



0. 
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FIG. 2. Same as Fig. [TJ but for the one-particle terms in the effective Hamiltonian V . 



This equation contains two different contributions, the first with total momentum P + + e 
and the second with total momentum P + . Projection onto a^(e)a^(P + ) |0) gives 



A 



1 XA e c 



2 



0. 



(2.27) 



and onto a)(xP + )a)({l - x)P+)|0) yields 



fi 2 + I + 



v 7 ! 



7T6 



x(l — X] 



A 



1 + -V(x) - rtAx) = 0. 

2 ^x(l - x) 2 1 ; 2i(1-j;) u ; 

(2.28) 



From the first projection f )2.27p . we conclude that A = or a = 1/2 ancfl 

u= -% = -(2^ 2 + /)/A. (2.29) 

V47T 

The valence eigenvalue problem is then trivially solved as M 2 = —ft 2 , and the auxiliary 
equation ( I2.28P reduces to 



- fi 



x(l — x) 



1 



A 



1 



t 1 (x) = 0. (2.30) 



Because the physical one-particle eigenstate is computed to have a negative mass squared 
of —fi 2 for any coupling value, one can construct multi-particle eigenstates of well-separated 
particles to have a Hamiltonian eigenvalue of any large multiple of — fi 2 . This is consistent 
with the demonstration by Baym [20] that the spectrum of <f) 3 is unbounded from below. 
Without the inclusion of zero-mode contributions, this would not be apparent in the LFCC 
method. 



1 We obtain the same result from a vacuum valence state and the action of (1 — P v yP on the vacuum. 
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III. 4> 4 THEORY 

The Lagrangian and light-front Hamiltonian density for 4 theory are 

1 / ^ , \ 9 1 9,0 A 



and 



£ = ^?--/<?- 1[ / (3-1) 



H = \^ + (3.2) 



The mode expansion for the field <ft is the same as (12 .3p for <ft 3 theory. We again split the 
light-front Hamiltonian V~ into two parts 

P free = / ^ + ^ (P + Hp + ) (3.3) 



and 



V ^ = ^S dx J^S 6{1 " ^ / hVM^>(* 2 p>(^ + ) (3.4) 

) a* (x2P + ) (X3P 4 " ) a (p + ) ] 



+ — / — — / a'(xp W((l — x)p )a(x p )a((l — x)p ). 

4 J y/ x (l - x)x'(l - x>) J 4tt v f ' vv " / Vv ^ ; 

A graphical representation is shown in Fig. [3j 



(b) 



FIG. 3. Same as Fig. [H but for </> 4 theory. Only the zero-mode term is included in (b). 

We focus on the zero-mode contributions to a vacuum valence state and consider the T 
operator 

r = i/>/i7rea t (e). (3.5) 



As before in (]2.10p . this corresponds to a shift of the field by a constant v. The relevant 



commutators with V are 

[Pta.ro] = *y^/iV(e), (3.6) 

+a) (xp + )a) ((1 — a;)p + )a(p + )] , 
[[^t,ro],T ] = ^ 2 A | ^aV)a(p+), (3.8) 

[[[Pr t , T ], T ], T ] = ^X^a\e), (3.9) 

where we keep only terms that will contribute to a projection 1 — P v onto states with one 
particle more than in the valence state. From these commutators we construct the effective 
Hamiltonian for the zero modes 

TMzm = u(jj? + u 2 X/Q)^a\e). (3.10) 
Figure H] provides a graphical representation. The first graph corresponds to the commutator 




- — -0 -^--© 




FIG. 4. Same as Fig. [21 but for t/> 4 theory. Only zero-mode terms are included.. 

in (I3.6P and the first term of V~\zm] the second to the triple commutator in (13 .9p and the 
second term of V~ |zm- Each commutator with V~ contracts one zero-mode creation operator 
with one annihilation operator in V~ . 

For a vacuum valence state, the valence eigenvalue problem is trivial, because P v V~\0) = 
0. The auxiliary equation fll.4p is projected onto the one-particle sector and becomes 

isi/j 2 + v 2 \/Q)J— a\e)\0) =0. (3.11) 



This implies that either v is zero or 



If we now consider the wrong-sign case, with jj? — > — /i 2 , we find v = ±v6A//i, which 
corresponds to the shift of the field <p that brings the Hamiltonian density to a minimum. 



fi 2 + \\v 2 = 0. (3.12) 
o 
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Thus, the inclusion of a zero mode in the LFCC T operator allows for the necessary shift 
in the field. Also, as can be seen from (13. 7p . the effective Hamiltonian will have terms that 
mix Fock states with odd and even numbers of particles, which is characteristic of broken 
symmetry. 



IV. WICK CUTKOSKY MODEL 



To illustrate what happens in a more complicated theory, we consider the Wick-Cutkosky 
model 19( of a charged scalar coupled to a neutral scalar. The Lagrangian is 



£ = \{d^f - l»V + (d, X f ~ m\ 2 - g<P\ X \ 2 , 



(4.1) 



where <fr is the neutral scalar field and \ the complex charged scalar field. The Hamiltonian 
density is 

U = ^ 2 + m 2 \ X \ 2 + g( p\x\ 2 . (4.2) 

The mode expansion for the field is the same as (12. 3 p for 3 theory; the mode expansion 
for x is 



X 



c + {p + )e 



(4.3) 



with c± the creation operator for the positive (negative) charge. The nonzero commutation 
relation is 

[c±(p + ),ci(p' + )]=5(p + -p' + ). (4.4) 
The free and interacting parts of the light-front Hamiltonian V~ are 



free 



dp^ 

p+ 



ji 2 a) \p + )a(p + ) + m 2 c^ + (p + )c + {p + ) + m 2 cL(p + )c^(p + ) 



2 J 



(4.5) 



and 



^int = 9 



dx 



'x(l — X) 



dp^ 



y/ATTp^ 



[a^(p + )c + (xp + )c_((l — x)p + ) 



(4.6) 



+ c+ (xp + )cL((l — x)p + )a(p + ) 
+ a\xp + ) \c\_((l — x)p + )c + {p + ) + cl((l — x)p + )c-(p + ) 

+ [c\(p + )c + ((\ — x)p + ) + cL(p + )c_((l — x)p + ) ) a(xp + ) 



A graphical representation is given in Fig. [5] 

We again focus on the zero-mode contributions and consider the T operator 

r „ = ^Si..(e) + W/dx* rf ,) e '(«)«L((l-«W. 



(4.7) 



The first term creates a neutral-scalar zero mode in the e — > limit; the second creates a 
neutral pair of charged zero modes. The first term again corresponds to a shift of the field 
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(a) 



(b) 



FIG. 5. Same as Fig.[3l but for the Wick-Cutkosky model. The charged scalars are represented by 
solid lines, with the arrow to the left (right) for positive (negative) charge; the neutral scalars are 
represented by dashed lines. The zero-mode contributions are labeled by v for the neutral scalar 
and by for the charged pair. 

by a constant v. The relevant commutators with V~ are 

[Vf Iee ,T ] = v\^a\e) + J ^^ )C \{xe)cl{{l - x)e\ (4.8) 



e 



[Pr t ,T ] = vg I f X cl(xe)cl((l-x)e)+g 



47re /3 f dxt (x) + . , 
-.a 1 e 



ev47re J yfx(l — x) 



+vg J ^( C t + (p+) c+(j9 +) + cUp + )c-(p + )), (4.9) 

[[Vto,To],T ] = 4*ug£ J ^^ c t.(xc) C l((l - x)e). (4.10) 

These yield an effective Hamiltonian of 

— — , /47r / 9 [ dxtn(x) \ +/ . , A 

p-i ZM ^_^ +9 _y_u_j » (E) ( 4 .n) 
+/* (*4£^<*> + ^r=y) c + (ie)c ' ((1 - I)£) ' 

Figure [6] shows a graphical representation. 

To have enough equations to solve for the unknowns v and t in the T operator, we 
must project the auxiliary equation (11.41) onto one and two-particle sectors orthogonal to 
the vacuum valence state. This yields 

e" f dxtn(x) 
vtf+g- —=2k=L= = (4.12) 



e 



y/x(l - x) 



11 



> 




FIG. 6. Same as Fig. El but for the zero-mode terms in the effective Hamiltonian V . 



and 



e P m 2 _|_ V g 



t (x) + 



0. 



(4.13) 



e x{l — x) 



x(l — x) 



A consistent solution is obtained if either v = and to = 0, which is a trivial solution 
without zero modes, or 



with (3 = 1, which shifts the field. We can speculate that this solution will lead to the 
expected unbounded spectrum. A complete calculation of the eigenstates would require 
addition of non-zero-mode terms to the T operator and calculation of various commutators. 
Because of the relative complexity of the Wick-Cutkosky model, compared to 3 theory, 
this is nontrivial and beyond the scope of the present work. However, we have shown that 
nontrivial zero-mode contributions to the T operator and, therefore, to the eigenstates can 
be computed in the LFCC approach. 

V. SUMMARY 

We have considered various examples of two-dimensional scalar theories where zero modes 
can play a role in the calculation of light-front Hamiltonian eigenstates. In each case, the 
LFCC method is able to incorporate the zero modes in a sensible way. To do this, the T 
operator must include terms that allow creation of modes with infinitesimal longitudinal 
momentum e, with the limit of e — > taken at the end. 

Although the examples are limited to two dimensions, there is no particular restriction 
on a direct extension to three or four dimensions. The zero-mode terms would include a 
dependence on transverse momenta. With these tools in place, one can use the LFCC to 
explore symmetry breaking nonperturbatively. 
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Appendix A: Singular terms 



For a more careful treatment of singular terms in the limit of e — >■ for <f) 3 theory, we 
begin with a more general zero-mode operator of the form 



T = J dp + g(p + )a\p + ), 



(Al) 



where g is limited to have most, if not all, of its support near zero. The commutators with 
V~ are 



[P^To] = J * + ^(p + )a T (p + ) 
[[Pr t ,T ],T ] = A J dp 



g{p + ) dp H 



y/4:np + J ^Jp' + (p + 

p+ g{p + ) dp'+ 



a' (p + )a(p + — p + ) 
g(p + -p' + )a\p + ), 



V^P + Jo \/p' + {p + -p' + ) 
[[Vr t ,To},To}\zu = J dp + ^g{p + )a\p + ), 

where I is the self-energy integral (I2.19p . The first term in (12.24ft becomes 

/i 2 + 1/2 A/2 



(A2) 
(A3) 

(A4) 
(A5) 



v- 



ZM 



dp A 



dp ,+g(p' + )g(p + -p H 



a\p + ^ 



p + V^P+ Jo ' ^p' + {p + -p' + ) 
The auxiliary equation (ll.4p . relative to a vacuum valence state, reduces to 

^ + 1/2 , A/2 [ p+ dpl+ g(P ,+ )9(P + -p' + ) =n 



(A6) 



p^ y/4np+ Jo ' ^/p' + (p+ -p' + ) 

We replace g by g{p + ) = g(p + )/ \fp~^ to obtain 



(A7) 



(/i 2 + -J)^ 



A/2 



dp' + g(p' + )g(p + -p' + ) 



(A8) 



To solve this equation, we introduce the Laplace transform G(s) = J °° e~ sp g(p + )dp + 
and make use of the fact that the Laplace transform of the convolution on the right-hand 
side is just the product of the transforms. This leaves 



(fx 2 + h)G(s) 



for which the solution is either G = or G 



A/2 ^ 



G 2 (s), 



(A9) 



^^■(li 2 + -/), a constant. The inverse 



Laplace transform of a constant is, of course, a Dirac delta function; therefore, we obtain 

9(P + ) = -4t?V+o WP + )> ( A1 °) 
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which, given the e — > limit to be taken, we can interpret as fixing g to be 



9(P + ) = "^V + - e). (All) 

This corresponds to the simpler expression ( \2.9\\ used for To in Sec. HIl with a = 1/2 and Aq 
given by the solution in (12.291) . 
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